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ABSTRACT 



GRB941017, a gamma-ray burst of exceptional fluence, has recently been 
shown to have a high-energy component which is not consistent with the stan- 
dard fireball phenomenology. If this component is the result of photomeson 
interactions in the burst fireball, it provides new and compelling support for sub- 
stantial high-energy neutrino fluxes from this and similar sources. In this letter, 
we consider what impact this new information has on the neutrino spectra of 
gamma-ray bursts and discuss how this new evidence impacts the prospects for 
detection of such events in next generation neutrino telescopes. 



1. Introduction 

Gamma-Ray Bursts (GRBs) are the most powerful objects in the universe, typically 
emitting luminosities between 10 50 and 10 54 erg/s over 0.1-100 seconds. Their 
energetics suggest that they may be the sources of the highest energy cosmic rays 
(Waxman 1995; Vietri 1995). It has been pointed out that, if high-energy cosmic rays are 
accelerated in GRBs, photomeson interactions between accelerated protons and target 
photons are inevitable, producing very high energy neutrinos and gamma-rays (see, for 
instance, Waxman & Bahcall 1997). Independent of any association of cosmic rays with 
GRBs, observable neutrino rates are predicted in models where similar energy goes into the 
acceleration of electrons and protons in the expanding fireball. 

In October of 1994, the Burst And Transient Source Experiment (BATSE) and the 
Energetic Gamma-Ray Experiment Telescope (EGRET) each observed an exceptionally 
powerful example of such an object. GRB941017 has the eleventh highest fluence observed 
in the nine years of BATSE observations. More interesting, however, is the fact that this 
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burst displays not only the typical synchotron-inverse Compton spectrum at ~30-1000 keV, 
but also shows a power-law high-energy component extending at least to 200 MeV in energy 
(Gonzalez et al. 2003). Other than this burst, EGRET has seen four GRBs at energies 
of ~ 100 MeV. These were each consistent with an extension of the synchotron-inverse 
Compton spectrum (Dingus 2001). Additionally, Milagrito has also observed evidence 
for emission near ~ 100 GeV in one burst. This observation lacked the ability to reveal 
significant spectral information, however (Atkins et al. 2000). In a burst-by-burst analysis 
of the complete BATSE catalogue, Guetta et al. (2003) identified GRB941017 as the most 
powerful neutrino emitter with, unfortunately, no neutrino telescope to observe it in 1994. 
Although other explanations may be possible, GRB941017 provides the best evidence to 
date for high-energy proton interactions with source photons in GRBs. 

If protons are accelerated to energies sufficient to produce the gamma-ray 
spectrum observed in GRB941017, high-energy neutrinos are a necessary consequence 
(Waxman & Bahcall 1997). Although no high-energy neutrino telescope of sufficient volume 
was operational in October 1994, it is interesting to consider the prospects for detection of 
a burst with similar characteristics in future experiments. 

In this letter, we calculate the neutrino spectrum predicted for such an event, and 
discuss the prospects for detection in future experiments such as the kilometer scale 
neutrino telescope IceCube (Ahrens et al. 2003). Our conclusions are: i) An event like this 
is likely to be observable with ~ 0.5 — 5 events from a single burst and, ii) a handful of 
events of this type can produce the diffuse flux of order 10 events per year predicted by 
fireball phenomenology (Guetta et al. 2003). 
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2. Neutrinos Prom Photomeson Interactions In GRB Fireballs 

Accelerated protons in GRB fireballs produce parent pions via the processes 

jry — > A — > n7r + (1) 

and, 

P7 -> A -> pvr (2) 

which have very large cross sections of a a ~ 5 x 10~ 28 cm 2 . The charged 7r's subsequently 
decay producing charged leptons and neutrinos, while the neutral 7r's decay into high-energy 
photons. To have sufficient center-of-mass energy for these processes to take place, protons 
must meet the threshold condition: 

4 >- (3) 

Primed and unprimed quantities refer to values measured in the comoving and observer's 
frames, respectively. In the observer's frame, 

r 2 

e p > 1.4 x 10 16 ^^eV, (4) 

^7,MeV 

resulting in a neutrino energy 

1 r 2 

e v = -(x p ^ )e p >7x l0 14 ^^eV, (5) 

4 £ 7 ,MeV 

where T 2 .5 = T/10 2 5 is the bulk Lorentz factor and £ 7i m c v — £ -y/ lMeV is the typical target 
photon energy. (x p ^ n ) ~ 0.2 is the average fraction of energy transferred from the initial 
proton to the produced pion. The factor of 1/4 comes from the assumption that the 4 final 
state leptons in the decay chain n + — > u^jjL + — > u^e + u e u^ equally share the pion energy. 

Typical GRBs display a broken power-law spectrum consistent with synchotron and 
inverse Compton emission. This spectrum can be parameterized by: 

e" a \ e y < e b 

F 7 = e 7 <in 7 /cfe 7 oc ^ 7 . (6) 
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In the case of GRB941017, a and f3 appear to change with time varying between 1.2 and 
2.1 for P and between —0.2 and 0.5 for a. In the first tens of seconds of the burst, the 
break energy is about 500 keV. It steadily drops, falling well below 100 keV after about 100 
seconds. 

For each proton energy, the resulting neutrino spectrum traces the GRB photon 
spectrum, but with a much higher energy break: 

4 = 7x10" 1 J 15 eV. (7) 

V ' ) fc 7,MeV 

z is the redshift of the GRB. For this particular burst, the redshift has been estimated to 
be z ~ 0.6, although its precise value is uncertain (Guetta et al. 2003). 

At very high energies, pions lose energy by synchrotron emission before decaying. This 
affects the spectrum of neutrinos from pion decay above the energy: 

< M = L~% 2 rUv,-2 eV, (8) 

where L 7)52 is the gamma-ray luminosity of the burst in units of 10 52 erg/s. For GRB941017, 
L 7)52 ~ I. t v -2 is the time scale of fluctuations in the GRB lightcurve in units of 10~ 2 s 
Neutrinos from muon decay have an energy cutoff that is 10 times smaller due to the longer 
lifetime of the muon compared to that of the pion. Above this energy, the slope of the 
neutrino spectrum steepens by two. 

The fraction of energy in accelerated protons converted to pions in the fireball is 
estimated from the ratio of the size of the shock, AR', and the proton mean free path: 

fir — ~ 7 (pp— »7r) ■ (9) 
Ap 7 

Here, the proton mean free path is given by A P7 = 1/ti 7 <7a, where n 7 is the number density 
of photons. The photon number density is given by the ratio of the photon energy density 
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and the photon energy in the comoving frame: 



n 7 



£' 



Using these equations, and the relationship i? ~ 2T 2 ct v , it is found 
which leads to 

Thus far, this calculation describes protons at the break energy. In general, 



£ I \ no -^7,52 j ( £ p/ £ p) a £p>£ p 

f n (e p ) ~ 0.2 7 ' x <^ , (13) 

1 2.5 r «,-2£ 7) MeV [ {e p /e b p f Bp < e\ 

where e b p is given by Eq.(4). 

This equation indicates that f n and, therefore, the spectrum's normalization and 
the final event rates, vary significantly from burst-to-burst (Halzen & Hooper 1999; 
Alvarez- Muhiz et al. 2000). Such fluctuations are constrained, however, and similar 
conclusions are reached when fixing /„. = 0.2, its typical value (Guetta et al. 2001). 

To estimate the bulk Lorentz factor, T, we relate it to the peak energy of the gamma-ray 
spectrum: 

£-1/2 

4 » r2 7 f MeV - ( 14 ) 

1 2.5^,-2 

For the case of GRB941017, we estimate T ~ 250 and this is again somewhat uncertain 
(Guetta et al. 2003). 

Finally, we obtain the neutrino spectrum: 

v de v ~ 81n(10) J " 1 ] 
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where, F 7 and f w are determined by equations 6 and 13, respectively. The neutrino 
spectrum found for GRB941017 in the time bins in (Gonzalez et al. 2003), is shown in 
figure 1. 

3. Event Rates in Neutrino Telescopes 

Currently, the AMANDA-II detector (Andres et al. 2001), located at the 
South Pole, provides the strongest limits for high-energy neutrinos from GRBs 
(Barouch & Hartdke 2001; Stamatikos et al. 2003). Considerably larger experiments such 
as the kilometer-scale neutrino observatory IceCube (Ahrens et al. 2003), are presently 
under construction. For a review of high-energy neutrino astronomy, see Halzen & Hooper 
2002 or Learned k Mannheim 2000. 

High-energy neutrinos are observed as muons or showers created in interactions inside 
or near a Cherenkov detector embedded in an optically transparent medium such as ice 
or water. Muons, created in charged current interactions travel several kilometers before 
losing the majority of their energy, thus enhancing their prospects for observation. The 
probability for observing an energetic muon is given by 

P v ^ = NaR^cc, (16) 

where is Avogadro's number, is the muon range (Dutta et al. 2001) and a cc is the 
neutrino-nucleon charged current cross section. For GRB941017, with a zenith angle of 
94.5° (for experiments at the South Pole), there is sufficient ice to allow for very long muon 
ranges. In our calculations, we take into account the absorption of neutrinos in the Earth 
as well as the effect of oscillations. 

All three flavors of neutrinos may interact producing showers within the detector 
volume. The probability of making such an observation is similar to the relation for muons, 
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although shower events do not benefit from long muon ranges, and typically they have to 
be produced inside the detector to be observed. More details on how the rate is calculated 
can be found, for instance, in Appendix C of Guetta et al. 2003. 

The rates for a GRB similar to GRB941017 are shown in table 1. We have made two 
estimates of the rate. The first estimate uses Eq. 15 for the neutrino flux, obtaining F 7 from 
a broken power-law fit to the observed gamma-ray spectrum. It is important to remark 
that this fit cannot accommodate the late-time high-energy feature of the gamma-ray 
spectrum. We predict 0.3 events per square kilometer, or roughly 1 event in IceCube 
with an effective area for muons and showers exceeding 1km 2 . Bear in mind that for an 
observation over a narrow temporal and angular window in coincidence with the optical 
display, the background is entirely negligible (Guetta et al. 2003). A single high-energy 
neutrino event is a conclusive observation. This rate is, of course, subject to the ambiguities 
associated with z and V. 

To accommodate the late-time, high-energy feature seen in the gamma-ray spectrum, 
instead of using the method described above, we can, alternatively, estimate the neutrino flux 
and event rate by relating the observed high-energy gamma-ray component of GRB941017 
to a corresponding high-energy neutrino component. This method assumes that the 
observed high-energy gamma-ray component is the result of photomeson interactions. 
The observed high-energy gamma-ray component of GRB941017 can be parameterized 
(Gonzalez et al. 2003) by 



dK 



7 ^ A 7 E~ S (17) 

where 5 is measured to be 0.96 to 1.10 (Gonzalez et al. 2003) and A 7 is a normalization 
constant fixed to the observed spectrum. The corresponding neutrino flux can be calculated 
by energy conservation given that the physics of photomeson interactions fixes the 
ratio of energy going to photons and neutrinos ( Alvarez- Muniz & Halzen 2003). The 
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main uncertainties are the slope of the neutrino spectrum, and the energy to which the 
high-energy gamma-ray component extends. Without cascading of the gamma rays in the 
source both would be uniquely determined by the photon spectrum. This is not to be 
expected and we therefore calculate the neutrino flux varying the neutrino spectral slope 
and the high-energy cutoff of the gamma ray spectrum. Note that the latter is unobservable 
because high-energy gamma rays are absorbed on diffuse infrared background radiation. 

The results are shown in figures 2 and 3. For a E' 1 neutrino spectrum and a 
gamma-ray cutoff near 100 TeV, we predict ~ 10 events per square kilometer. Experiments 
capable of observing TeV gamma-rays from GRBs, such as MILAGRO (Atkins et al. 2001) 
and IceCube (Halzen & Hooper 2003) will be useful in determining the maximum energy to 
which the photomeson component of the gamma-ray spectrum extends. 

Another way of summarizing our results is shown in figure 3 where the rate is plotted 
as a function of the neutrino spectral slope. The observed event rate increases with the 
spectral index because more neutrinos of lower energy are produced and, although their 
detection probability is smaller, their larger number compensates. This trend continues 
until a neutrino spectral slope of ~ 2.1. Above this slope, a significant amount of the energy 
goes into neutrinos below the experimental energy threshold of the detector and the rate 
decreases. The same argument explains why for a flatter neutrino spectrum, similar to the 
observed photon spectrum, the same event rate is reached for a higher gamma-ray cutoff as 
seen in figure 2. 

Although ~ 10 3 GRBs are observed per year, the total rate from all of these events is 
typically estimated to be ~ 10 events per square kilometer per year (Guetta et al. 2003). 
A small number of exceptional bursts, such as the one we study here, can contribute 
substantially to the total neutrino flux. 
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4. Conclusions 

We have pointed out that the high-energy feature observed in GRB941017 provides 
further support for the expectation of detectable fluxes of high-energy neutrinos in 
coincidence with gamma ray bursts. We have estimated the neutrino spectrum that would 
accompany such a burst and discussed the prospects for such an event's detection in future 
neutrino observatories. We find that such an event is expected to produce on the order of 1 
event in a kilometer scale neutrino telescope and that this would be a conclusive observation 
since there is no competing background during the time and in the direction of the burst. 

J. Alvarez-Muhiz is supported in part by Xunta de Galicia (PGIDT02 PXIC 20611PN) 
and by MCyT (FPA 2001-3837 and FPA 2002-01161). F. Halzen work is supported in part 
by DOE grant No. DE-FG02-95ER40896, NSF grant No. OPP-0236449 and in part by the 
Wisconsin Alumni Research Foundation. D. Hooper is supported by the Leverhulme Trust. 



-Il- 



ea 
o 



> 



> 

O 

1 

CM 

W 



10" 



10" 



10" 



10" 



10 



10 



10 



12 



TTTTT 1 1 I I I I III 1 1 I I I I III 1 1 I I I I III 1 1 I I I llll 1 1 I I Mill 1 1 I I I I III 1 1 I I I I II 



X 



* x □ 

jr □ 

□ 

□ 

a ■ 

□ ■ 
□ ■ 

.□ ■ 



X 

X 

x* x~ □ 
x x x □ 
x x A □ 
**x n 
*? x D 

□ _ 



x- x S Rxxxxx^ xxxxxxx V 

X X m a ■ X *[ 



[] 



X 
X 

□ 



X *fl B 

x mm 

x *□■ 

X SK 1 



X ^ n 

x *□ ' a 
x x □ " B 

x x □ I 

x x □ 
X x Q 

X X 



□ 

□ 



x X p 

X X □ 1 
X * □ 
X * □ 
X * □ 
x * $ 
X 

X *J 

X ! ! 

X 

X 



■ Time bin 1 
Time bin 2 
Time bin 3 
Time bin 4 ■ 

imJ i I i I i I i I i I i I i 



10 3 10 4 10 5 



10 6 10 7 
E v [GeV] 



10 c 



10 9 10 10 



Fig. 1. — The high-energy neutrino spectrum estimated for GRB941017. The spectrum is 
shown for the four time bins in Gonzalez et al. 2003, the first of which begins 18 seconds 
prior to the BATSE trigger time. Each bin is of ~33 second duration. In Gonzalez et al. 
2003, a fifth time bin is described which we find to have a negligible impact on the neutrino 
spectrum. 
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Maximum Energy of y of hadronic origin [TeV] 

Fig. 2. — The muon neutrino event rate per square kilometer in a high-energy neutrino 
telescope for GRB941017 as a function of the energy to which the high-energy gamma- 
ray component observed by EGRET extends. The rate is calculated relating the observed 
high-energy gamma-ray component of GRB941017 to a corresponding high-energy neutrino 
component. Results are shown for two choices of accelerated neutrino spectra: E~ x (solid 
line), E~ 2 (dashed line). Protons are assumed to be accelerated up to E = 10 20 eV. The 
muon energy threshold is E* hr =100 GeV. 



-13- 



> 

e 
o 
3 



10 



1 - 



0.1 - 



0.01 



- 1 1 1 1 ! 

E^OO GeV 
E^lOO GeV 



1.5 2 2.5 

Neutrino spectral index y, dN/dE=A E" Y 



Fig. 3. — The muon neutrino event rate per square kilometer in a high-energy neutrino 
telescope for GRB941017 calculated by relating the total energy of gamma-rays in the high- 
energy component to the total energy in high-energy neutrinos. Results assume the high- 
energy gamma-ray component observed by EGRET extends to 1 TeV. Protons are assumed 
to be accelerated up to E = 10 20 eV. The effect of energy threshold is also shown: E* hr =500 
GeV (solid line); £* hr =100 GeV (dashed line). 
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Calculation Method 


Muons 


Showers 


Calculated From Target Density 


/tt Calculated 


0.27 


0.039 


fn Fixed (=0.2) 


0.091 


0.013 


Inferred From Energy Conservation 


^max = 10Q GeV 


0.21 


0.065 


£max = X Te y 


2.1 


0.65 


£max = 1Q Te y 


21. 


6.5 



Table 1: The event rate (per square kilometer) from GRB941017 for two calculational meth- 
ods. Those rates found by the target density calculation (as described in Section 2) are 
given for the calculated value of f n the fraction of proton energy converted to pions, and for 
a fixed value of /„- = 0.2. This calculation assumes that roughly equal amounts of energy 
go into accelerating protons and electrons. The rates inferred from energy conservation are 
shown for three different choices of the maximum energy to which the high-energy gamma- 
ray component of the spectrum potentially extends. An E~ 2 proton injection spectrum is 
used. Both muon and shower energy thresholds are set to 100 GeV. 
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